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Abstract 

We have calculated the free energy up to two loop to compare T 2 with T 4 in IIB matrix model. It 
turns out that T 2 has smaller free energy than T 4 . We have also discussed the generation of the 
gauge group by considering fc-coincident fuzzy tori and found that in this case U(l) gauge group 
is favored. This means that if the true vacuum is four-dimensional, it is not a simple fuzzy space 
considered here. 
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1 Introduction 



IIB matrix model £Q is a candidate of the constructive definition of superstring theory. It is defined 
by the action 



where and \& are N x N Hermitian matrices. is a ten-dimensional vector and \& is a ten- 
dimensional Major ana- Weyl fermion. 

One of the important aspects of this model is the dynamical generation of the spacetime and the 
gauge group. In this model the eigenvalue distribution represents spacetime geometry. From this 
point of view, if a four-dimensional eigenvalue distribution has smaller free energy than the other 
configurations, the spacetime will be compactified to four-dimensions. On the other hand the stability 
of k coincident D-branes may indicate the dynamical generation of U (k) gauge group. So far, various 
attempts on this have been made. Several such attempts using branched polymer [2], complex phase 
effect and especially improved mean-field approximation 0] were successful. In these studies, 
four-dimensional configurations were found to have smaller free energy. 

Another interesting approach is to compare the free energies on noncommutative backgrounds. In 
string theory, noncommutative gauge theories on flat or curved backgrounds are realized with con- 
stant or nonconstant 0. IIB matrix model on flat fuzzy space backgrounds can be mapped to 
noncommutative super Yang-Mills theory (NCSYM) 0, and the system can be analyzed by pertur- 
bative calculations. Noncommutative gauge theories on curved manifolds have recently been studied 
extensively by deforming IIB matrix model |S| . 

However, flat fuzzy space backgrounds cannot be realized as classical solutions at finite N. There- 
fore, in order to see the iV-dependences of the free energies, we must study compact noncommutative 
spaces. Fuzzy S 2 , S 2 x S 2 and S 2 x S 2 x S 2 backgrounds have already been studied and S 2 x S 2 has 
been shown to be stable under certain deformations E3j • 

In this paper, we study the free energies of IIB matrix model on fuzzy torus backgrounds. The 
relation to flat fuzzy space is more transparent in the case of the fuzzy torus than in the case of fuzzy 
sphere. We expand the action of IIB matrix model around a fuzzy torus background and sum up all 
1PI diagrams. Then, we search the minimum of the free energy. Especially, we compare T 2 and T 
backgrounds. In principle, in order to find the true vacuum, we must calculate free energies on all 
possible backgrounds and find the minimum. Here, however, we consider only fuzzy torus. Although 
fuzzy torus is not a classical solution, however, as we will show later, fuzzy torus locally looks like 
flat fuzzy space, which is a classical solution. Therefore, we expect our calculation captures some 
qualitative feature of flat fuzzy spaces. If the true vacuum is similar to a flat fuzzy space, we may find 
some informations such as dimensionality of the spacetime on it. In present study we will see that T 2 
has smaller free energy than T . This implies that if the true vacuum is four-dimensional, it is not a 
simple noncommutative space considered here. 

Although both fuzzy sphere and fuzzy torus break supersymmetry, the breaking effect is soft and 
the loop corrections to their free energies almost vanish. This is because their local structures are 
similar to flat fuzzy spaces. In the case of fuzzy S 2 , S 2 x S 2 and S 2 x S 2 x S 2 , such cancellations have 
been confirmed to 2-loop level Elj ■ For fuzzy torus we can see such cancellation at any loop level. 
Furthermore, in this case the calculations become simpler and we can study larger parameter region. 
For example, on T 4 , we vary the ratio r of the noncommutativity parameters of two tori and their 
ratio R of the matrix size 1 . We can also see the effect of the global topologies on the free energies by 
comparing the results for fuzzy tori and fuzzy spheres. 

On the other hand, fuzzy torus has a few drawbacks. We embed fuzzy T 2 into R 4 , while fuzzy S 2 
can be embedded into R 3 . Therefore, in ten dimensions we can realize only two- and four-dimensional 
configurations using fuzzy tori ( T 2 and T 4 ), while a six-dimensional configuration can be constructed 
using fuzzy spheres ( S 2 x S 2 x S 2 ). In ^0 it is shown that S 2 x S 2 has a smaller free energy than 
S 2 x S 2 x S 2 . In our case, such comparison is impossible. We also have subtlety in the perturbative 

1 In the case of S 2 X S 2 , r 2 = R was studied 9 . 




(1.1) 
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calculation, due to the presence of a few tachyons on fuzzy torus. We will discuss this problem in 
detail later. 

The organization of the paper is as follows. In section 2, we will construct the fuzzy torus back- 
grounds and show the correspondence between fuzzy torus and flat fuzzy space. In section 3, we 
calculate the free energy up to two loop level analytically, and in section 4, we compare the free en- 
ergies numerically. In appendices, we show the detail of the calculations. Although we are interested 
in T 4 , we will mostly show the Feynman diagrams and calculations of T 2 as T 4 calculations are very 
lengthy. 

2 IIB matrix model on fuzzy torus backgrounds 

In order to perform a perturbative calculation, we decompose and \& into a background and 
the fluctuation a M , ip around it: 



i„ , * = + (p. 
Adding a gauge fixing and the corresponding ghost term 



<Sgfh 



FP 



Tr -[P M ,a"] 2 + [P ll ,b\W 



(2.1) 



(2.2) 



to we obtain the full action 



4.9 2 



r Tr (V^a^F^ 
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'4^2 Tr [ 8 (" a J 2 ' 



(2.3) 



where Vy. = • ], F^ v = P„] and T^u = [F^u, ■ ]• If the background is not a classical solution, 
e.g. fuzzy sphere or fuzzy torus, the 0(a^) term does not vanish. 

2.1 IIB matrix model on fuzzy T 2 background 

First, we introduce the N x N clock and shift matrices U and V : 
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where U and V satisfy 
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Using them, we construct a 2-dimensional fuzzy torus as follows : 

Pi + iP 2 = 2\fNLtU, Pi - iP 2 = 2v r NL 1 U~ 1 , 
P 3 + iP A = 2VNL 2 V, P 3 - iPi = 2VNL 2 V~ x , 

p^Vnl^u + u- 1 ), p 2 = -iVNLi(u -u- 1 ), 

P 3 = VnL 2 (V + V- 1 ), P 4 = -iVNL 2 (V-V- 1 ), 

Ps = ■ ■ ■ = Pw = 0. (2.7) 

Here, L\ and L 2 are the noncommutativity parameters, which have the dimension 2 of mass. Their 
product, L 2 — L\L 2 , corresponds to the noncommutativity parameter B of fuzzy plane. The relation 
of this background to fuzzy plane is discussed more precisely in section \2. 41 

We can expand the fluctuations in terms of U and V. For example, a M can be expressed as 

N 

a„ = ^(m u m 2 )U m W m \ (2.8) 

mi ,T7l2- 1 

and the hermiticity condition is given by 

~a;(m u m 2 ) = e -2~/^ 5fi (_ mi) _ m2 ). ( 2 . 9 ) 

We will show that this expansion corresponds to the ordinary Fourier expansion in section 12.41 
Using these Fourier modes, we can construct the Feynman rules, which we summarize in Appendix 1X1 

2.2 IIB matrix model on fuzzy T 4 background 

We can similarly construct fuzzy T 4 using the tensor products of the fuzzy tori constructed above : 

P a = yfW x Li{U x + C/f 1 ) <8 1jv 2 , P 2 = -i^/WxLiiUx - C/f x ) ® 1 N „ 

P 3 = yfihLi(Vi + V,- 1 ) <g> 1 N2 , P 4 = -i^/NlLx{Vi - Vf 1 ) ® 1jv 2 , 
P 5 = V^L 2 l Nl ® (U 2 + U2 1 ), P 6 = -i^L 2 l Nl ® (Z7 2 - t/^ 1 ), 
P 7 = ^N~ 2 L 2 l Nl <g> (U 2 + U^ 1 ), P 8 = -iy/W 2 L 2 In, ® (U 2 - V^ 1 ) , 
P 9 =P 10 -0. (2.10) 

Here, A~ = NiN 2 , and E/j, V* are the A/j x TV^ clock and shift matrices. For simplicity, we set the radii 
of the two circles in the first T 2 to the same value 2*/N\L\, and similarly in the second T 2 to 2\fW 2 L 2 . 
In the following sections, we first fix the value of N, and calculate the free energy as a function of 
R = j^, L = \[L~\L 2 and r = j^, and minimize it with respect to them. 

2.3 /c-coincident fuzzy tori 

We also consider a system of fc-coincident fuzzy tori. To study fc-coincident fuzzy T 2 , we replace Urm 
and Vrm with U(N/k) ® lfc an d ^(JV/fc) ® Ifcj where Urm and V(jv) stand for the N x N clock and shift 
matrices, respectively. Then, the Fourier modes a^,(f,b and c become fc x fc matrices. 

Similarly, in the case of fuzzy T 4 , we can replace Ui ® lj\r 2 with [/1 (8 1jv 2 (8 lfc, and so on, to 
construct a fc-coincident T 4 . In this case, we have N = N\N 2 k. 

2 In IIB matrix model, usually we treat as the spacetime coordinate. However, when we map IIB matrix model 
to the non-commutative Yang-Mills theory, we usually regard P M as the momentum operator and assign the mass 
dimension to A^. 
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2.4 Correspondence between fuzzy torus and fiat fuzzy space 



In this section, we discuss the relation between the fuzzy torus and flat fuzzy space To be 

specific, we consider the case of fuzzy T 2 and fiat fuzzy plane ( flat Dl-brane ). Generalization to 
higher dimensions is straightforward. 

In order to avoid confusion, we denote fuzzy plane background as Q M and that of fuzzy torus as 
Pp. The fuzzy plane is characterized by the constant commutators : 

[QuQi] = -iB, Q 3 = ...=Q w =0. (2.11) 

Let us write U and V as 




U = exp^^Q[\, V^exp[ lx jif E Q' 2 ). (2.12) 

Then, if we consider the region where eigenvalues of Q' ± and Q' 2 are close to zero, the condition 
UV = e 2wi / N VU is equivalent to [Q[,Q' 2 ] = -iB. Therefore, Qi ~ Q[ in this region. On the other 
hand, if we compare 12. 7|) with 



11 ~ 2 



yfW>. ,, 13) 



we have the identifications in the embedded space as P2 — const. Q'^P^ ~ const.Q' 2 . In general, 
tangent spaces of fuzzy torus look like fuzzy planes. 

IIB matrix model on fuzzy torus can be mapped to NCSYM as in the case of flat fuzzy space 
[HI E] • If we introduce the spacetime coordinates as 

xi = --g-, x 2 = -g, (2.14) 



we have 



= e^*^* X (phase), (2.15) 

where momenta fcj are defined by fcj = rrii^j ^jp-- Hence the expansion J23 becomes the ordinary 
Fourier expansion in the commutative limit. 



3 Calculation of free energy 

We calculate the free energies in this section. We ignore the tachyons and zero modes on the fuzzy torus 
in our calculations. The vector- and fermion-kinetic terms on fuzzy torus have tachyons and extra 
zero-modes besides those coming from the U(l)-part. (See Appendix iDl for the detail.) Therefore, 
it is not clear whether perturbative calculation is possible or not. However, the number of tachyons 
remains the same in the large- N limit, and their mass squared is of O (-4) and vanishes in the large- TV 

limit. Hence we expect that the tachyons give a relative correction of O ( '"^^ J to the free energy. 

On the other hand, in the case of T 4 and Ni ~ JV2 ~ V~N, we will see that the leading contributions 
cancel out due to the supersymmetry, but the terms of O (-^-) do not vanish. They are expected to 
be larger than the effect of the tachyons, and we can take the large- N limit without paying much 
attention to the tachyons. 
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3.1 Tree-level free energy 



We start with the tree-level free energy. It is given by 



1,3 



^ 2 ~ sm (n) 

16Ntt 2 L 2 



for T 2 , and 



16L 2 iV 3 , / 7T 



/ 7T N 

(in the large- A limit) 



(3.1) 



■ sm 



16Nn 2 L 2 
9 2 

for T 4 , respectively. 

3.2 1-loop free energy 

The 1-loop free energy can be easily evaluated as pQ 



r 2 H — - ) (in the large- N limit) 



(3.2) 



2 ^ n 

n—1 



1 / 2i \ n 1 ^ 1 / i 



4 ^ n 

n— 1 



y7' 



(3.3) 



The terms for n — 1,2,3 cancel out, and the first nonzero contribution comes from n = 4. In the case 
of T 4 and Ni ~ ~ \A/V, it is 0(N°) and can be neglected compared with the tree- level and 2-loop 
level contributions. In the case of T 2 , truncation at some order of T is not a good approximation 
because of infrared divergence. 

Instead, we can numerically evaluate the 1-loop free energy without expanding in T . We calculate 
all the eigenvalues of the kinetic terms, and sum up their logarithms. We find that the 1-loop free 
energy on T 2 background is proportional to log N and hence negligible compared to the tree level and 
2-loop contributions. 3 



3.3 2-loop free energy 

Next we evaluate the 2-loop free energy on the fuzzy T 4 . Since it is difficult to obtain analytic forms 
of the vector and fermion propagators, we expand them in powers of T as shown in Appendix IA.1I 
and calculate the free energy order by order in T . A naive power counting shows that each power of 
T is associated with and ^= for T 2 and T 4 , respectively, if there is no infrared divergence. In fact, 

this expansion is good on the T 4 background, while it is not on the T 2 background. We discuss the 
validity of this expansion in Appendix [C] Here we show the result for the T 4 background. 
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3.3.1 0(T°) 

We have four kinds of 1PI 2-loop vacuum diagrams ( see FigQJ, and their contributions can be 
approximated as (see Appendix [E] for the detail) 

g 2 V l-cos(27rmon) 
[ ' 87V 2 L 4 ^ F(m) 2 -F(n) 2 ' 1 ; 

ff 2 y 1 ~ cos(2^m o n) 



512^ 2 L 4 ^ F(m) 2 • F(n) 2 

m,n v ' v ' 

n 27 g 2 \ - 1 - cos(27rm o n) 

[C) ^ 512N 2 ^ F(m) 2 -F(n) 2 ' 1 ' 



(A\ - 45 g \ - 1 - cqs(2ttto o n) 

1 ' ~* 256N 2 L> Flm) 2 ■ Fin) 2 ' 1 ' ' 



where mon = m^ v n v = miM^lam + m 3 n 4 -m 4 n 3 ^ &nd p{jl y = ^\ F ^ . 

(n) are analogous to the ordinary momentum variables in the continuum limit and given in appendix 

M 

Summing up these approximated values H3.4|) ^ (|3.7[l . we nave 

I__L 27 45 \ g 2 v 1 - cos(27rm o n) 

8 512 512 256 I N 2 ^ F(m) 2 ■ F(n) 2 { ' 

This cancellation is not exact, but we can see numerically that the exact contribution of O(J r0 ) 
terms is smaller than that of 0(!F 2 ) terms. 

3.3.2 OiT 1 ) 

The O^ 1 ) contribution is zero, because integrands are odd function of the momenta. 

3.3.3 0{T 2 ) 

The calculation of 0{T 2 ) contribution is very tedious. Here we show only the simplest case, that is, 
the planar part of Fig^ (b). In this case, the contribution can be writen as 

q 2 N 

9 —h{R,r), (3.9) 



3 In this calculation, we have discarded the tachyons and the extra zero-modes. 
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where 



fb(R,r) 



d 4 kd 4 l 



1 - cos (27r(fci + h)) cos (27r(fc 2 + Z 2 )) 



256 J { F(k) 2 ■ F(l) 2 ■ (F(k + l) 2 ) 3 

x = (sin(7rfci) sin(7rZi) cos (7r(fci + h)) + sm(nk 2 ) sin(7r? 2 ) cos (7r(fc 2 + h))) 

r 3 \/R 

1 - cos {2Tr(k 3 + l 3 )) cos (2n(k 4 + Z 4 )) 
+ F(k) 2 -F(l) 2 -(F{k + l) 2 ) 3 

xr 3 v / i? (sin(7rfc 3 ) sin(7rZ3) cos (7r(fc 3 + lg)) + sin(7rfc4) sin(7rZ 4 ) cos (7r(fc 4 + Z 4 ))) 

(3.10) 

Here A;, = m,i/Ni,kj — m^/JVa, lj = n,/iVi,Zj — rij/Nz (i = 1,2; j = 3,4), and fc M and Z M vary 
from to 1. Because the ghost propagators do not contain JF, this expression is rather simple. 
The contributions from the other diagrams ( Fig.l (a),(b),(c) ) also can be written as ^jj-f a {R, r), 

2jj-fc(R,r) and a rr-fd(R, r). But the expressions of f a (R,r), f c (R,r) and fd(R, r) become much 
more complicated, because both vector and fermion propagators contain 0(T) corrections. Summing 
up all of the diagrams, we find that the leading behavior of the 0{T 2 ) contribution can be expressed 
as 

g 2 N 



L 4 



h(R,r). (3.11) 



The function h(R, r) = Y]j— a & c d fj(R' r ) * s dimensionless and has a symmetry under the exchange 
of two tori 

h(R,r) =h(R-\r- 1 ). (3.12) 

3.4 Total Free energy 

The total free energy on the fuzzy T" 4 background can be written as 



where we have neglected subleading terms in N. The 0(g ) terms come from higher-loop diagrams. 
If we neglect them and minimize the free energy with respect to L, we obtain the minimum value 





E{R 1 r-N)=ATTN x \[r 2 + — ) h(R,r) (3.14) 
at 

r4 1 / hlu „\ 

(3.15) 

The situation does not change qualitatively, even if we take the full order of the loop expansion 

(2 \ H— 1 
j-z) N 2 . However, the O^ ) 

contributions cancel due to the supersymmetry 4 , and the (^(J- 1 ) contributions vanish because they 
are given by integrals of odd functions. Therefore, if we include the contribution from the higher-loop 
diagrams, the free energy can be expressed as 

L 4 



71 (3-16) 



where hj are of 0(1). The minimum value is again of O(N). 

4 See Appendix IbI for the detail. 
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3.5 The case of ^-coincident fuzzy tori 

Here we consider the system of fc-coincident fuzzy tori. This system can be mapped to V(k) NCSYM. 
The free energy in this case becomes 

kN Y, (it) h ^ r )- ( 3 - 17 ) 

j=0 ^ ' 

Therefore, the single fuzzy torus has the lowest free energy, and U(l) gauge group is favored in this 
case. 



4 Numerical comparison of free energies 

As the analytic expressions are very complicated, here we try to compare the free energies numerically. 
First, we consider h(R,r). It is given by an 8-fold integral, and we have evaluated it approximately 
using the quadrature by parts. 

We divide the integration domain into M bins in each direction, with M = 9, 10, • • • , 15. Let 
h(R, r; M) be the value of h(R, r) evaluated by M bins. We have fitted the plot of h(R, r; M) versus 
1/M by a polynomial, and estimated h(R, r) = h(R, r; oo) extrapolating it up to 1/M 0. In Fig[21 




0.05 0.06 0.07 0.08 0.09 0.1 0.11 0.12 

1/M 



Figure 2: h(R,r; M) versus 1/M for R = r = 1. The curve shows a nice polynomial fitting. 



we plot h(R,r;M) against 1/M for R = 1 and r — 1. The curve shows a nice polynomial fitting. We 
can evaluate the h(R, r) = h(R, r; oo) replacing 1/M — in the polynomial. 

Once we have h(R,r), we can readily obtain E(R,r;N) from eqn. I|3.14|l . We plot E(R = 
l,r;N)/N as a function of r for R fixed to 1 in Fig0 We see that R = r = 1 is the minimum. 
Therefore, the fuzzy T 4 is favored under the variation of the ratio r of the noncommutativity. 

FigQ] shows E(R, r = 1; N)/N as a function of R for r fixed to 1, and this time R = r = 1 is the 
maximum. Therefore, we find that the fuzzy T 4 is unstable under the variation of the ratio R of the 
matrix sizes. 

To compare with the result for the fuzzy S 2 x S 2 in [S], we plot E(R, r — \J~R; N)/N as a function 
of r in Fig|SJ The curve r = VR is the same as the one studied for the fuzzy S 12 x S 2 in 0. On this 
curve, R = r = 1 is the minimum, hence the fuzzy T 4 is favored. This observation is similar to the 
fuzzy S 2 x S 2 . In HJ, it is pointed out that although the 2-loop contribution (h(R,r — \/~R) in our 
notation) is maximum at R = r = 1, the tree-level contribution, ~r 2 + ^ = i?+-^,is minimum at 
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35 1 1 1 1 1 1 1 1 1 

0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

r=L 2 /L 1 



Figure 3: E(R,r; N)/N against r for fixed R = 1. i2 = r = 1 is the minimum here. T 4 is favored 
under variation of r. 



38.2 




R=N 2 /N 1 



Figure 4: E(R,r;N)/N , against i? for fixed r = 1. R = r = 1 is the maximum here. T 4 unstable 
under variation of R. 



R = r = 1 and this makes the total free energy minimum at R = r = 1. We have exactly the same 
situation here. 

38.5 




R=N 2 /N 1 



Figure 5: E(R, r; N)/N against R for R 



R 



1 is the minimum here. T is favored. 



These results show that although the fuzzy T 4 is favourable in some directions, it is unstable under 
the variation of R. Because the limit R — > or oo corresponds to the fuzzy T 2 , we conclude that the 
fuzzy T 2 is more favourable than the fuzzy T 4 . Since our observations along r = VR is completely 
analogous to that of [Hj , it is plausible that the fuzzy S 2 x S 2 also has the same kind of instability. A 
similar phenomenon is observed in a 6 dimensional bosonic matrix model with Chern Simons term, 
where S 2 is indeed shown to be more favourable than S 2 x S* 2 |12| . 

We can compare the free energy of the fuzzy T A and the fuzzy S 2 x S 2 . At R = r = 1, the free 
energy of T 4 is ~ 38.1-ZV. According to [5], the free energy of S 2 x S 2 with the same matrix size and 
noncommutativity is ~ 3.6A. Therefore, S 2 x S 2 is more favorable than T 4 . 

We have also computed the 2-loop free energy of the fuzzy T 2 up to 0(T 2 ) as a function of the 
ratio of radii, L±/L2- We found that it has the minimum at L± — Li and is proportional to N. 
However, because the T expansion is not good in this background, this result is not reliable. The free 
energy on the S 2 background is calculated in ^H]- In this case, the tree-level action dominates and 
the free energy is at L = 0. 



5 Conclusions and discussions 

In this paper, we have calculated the free energies of IIB matrix model on the fuzzy T 4 backgrounds 
up to 2-loop level. The free energy on T is a function of the ratio r of the noncommutativity and 
R of the matrix size. At R = r = 1, the free energy takes minimum when r is varied, but maximum 
when R is varied. This indicates that T 2 is more stable than T 4 . It is plausible that S 2 x S 2 has the 
same sort of instability. If so, a four-dimensional vacuum might not be realized as a direct product 
of two noncommutative surfaces. Commutative backgrounds seem to be more plausible. Another 
possibility is that new terms like Myers term arise dynamically (for example via large- TV RG ) and 
stabilize noncommutative spaces. 

We also have compared the free energies of S 2 x S 2 and T 4 , and found that the latter is about 
10 times as large as the former. It implies that IIB matrix model favors configurations with higher 
symmetry. 

Furthermore, we calculated the free energy of the system of /c-coincident fuzzy T , and found that 
the k = 1 system has the smallest free energy. In terms of NCSYM, this means that U(l) gauge group 



10 



is favored. This result is the same as that of the fuzzy sphere in IIB matrix model |5| and in bosonic 
matrix models with Chern Simons term |13j . It seems that in general noncommutative backgrounds 
favor U(l) gauge group. 
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A Feynman rules of IIB matrix model on fuzzy T 2 background 

In this section, we summarize the Feynman rules of IIB matrix model on the T 2 background. Gener- 
alization to the T background is straightforward. 

A.l Propagators 
Boson ( vector and scalar ) 

The kinetic term for is given by 

-^■Tr a M {V 2 S f _ ll/ - 2iT ilv )a v 



and the propagator is expressed as 



(A.l) 



M-mi.-ma^m.na)) - 9 - {VH^ - 2^}^.^ x e a*i»W". (A.2) 
We expand it with respect to T and truncate it at some order : 

Sftu . 1 1 1 1 1 1 

Note that this expansion is not good in this case but a similar expansion in the T 4 background is 
reliable. 

Fermion 

The fermion kinetic term is given by 

- ±Tr ^fV = -^£ e- 2 ™^/ N <p T (- mil -m 2 )T mim2 ., nin Mni,n 2 ), (A.4) 

and the propagator is expressed as 

(^•(-mi,-m 2 )^(nx,n 2 )) = ^(^) * e 2 ™^. (A.5) 

ij;mim2;nin2 



N \ f 
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Here we also expand i in powers of T 



-l 



(A.6) 



Ghost 

The Ghost kinetic term and propagator are given, respectively, by 

4jTr bV 2 c = % Y e- 2mm ^/ N b{~ mi , -m 2 ) (V 2 ) c(m,n 2 ), 

q2 q2 Z ✓ v ' V /mim2;ni?T,2 v y 



(A.7) 



and 



(&(-«!, -ma)c(n 1|n3 )) = | (P 2 )^;^ x e 2 — (A.8) 

A. 2 Explicit form of V 

A direct calculation shows that 

[PL,a„] = ^iVNLx^U^V-^ sin(^) {o P (mi,m2 - lje**™ 1 /* - a M (m 1 ,m 2 + i) e —Wivj , 
[P 2 , 0/t ] = £ C/^F" mi sin (^i) {o M (roi,ma - lje™" 1 ^ + a M (m 1; m 2 + l^"™"^} , 

m 

[P 3 ,aJ = -2 l V7Vi 2 ^C/ m2 F- mi sin(^) {a M (mi + 1, nuje"^ - a^rm - 1, mjje^} , 

m 

[P 4 , 0/ ,] = -2^iVL 2 ^?7 TO2 F- mi sin(^) [a^mi + 1, m^e"™" 12 ^ + a M (m x - 1, m^e™™ 2 ^} . 

m 

(A.9) 

From these equations, we obtain 

= f ±2z^i 1 sin(^) e ™ 1 /^ ((n 1 ,n 2 ) = (mi > m 2 ±l)) im 
Umim2 ' ni ™ 2 \ (otherwise) ' \ ■ > 

and so on. This expression is not very simple, but V 2 can be written in the following simple form : 
(r 2 ) mim = l^L 2 sin 2 (^) S mini S m2n2 . (A.ll) 

i 

A. 3 Interaction vertices 

There are four types of interaction vertices ( FigEl ) ■ Here, the solid, dotted and wavy lines represent 
the fermions, ghosts and bosons, respectively. 

(a) Fermion-fermion-boson coupling V3/ 
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(a) (b) (c) (d) 

Figure 6: Interaction vertices 



V 3f = <p i (k)('f) ij a li (-k - ?)^(?) e --{ fc i^+( 1 i2+(fc 1 +( 1 )(fc2+( 2 )}/iv sin(?rfc ^ 

g k,i 

Note that /x = 1, 2, • • • , 10 ( Both vectors and scalars couple to the fermion ). 
(b) Ghost-ghost-vector coupling V$ g 



(INY'i 1 

V 3g = 2i± 1 — J2e~™ {mim2+nin2+{mi+ni){m2+n2)}/N b(rn)c(n) 

^ m,n 

x sin (^i) e -^(m 1+ni )/N gin (^mon- ^) a 1+ {-m 1 - n u -m 2 - n 2 - 1) 
+Li sin (^i) e «(mi+m)/JV sin ^ m Q n + ™i ^ 5l _(_ mi _ ni) _ m2 - n 2 + 1) 
-L 2 sin e -^(m2+n 2 )/iV gin ^ mon+ ^ 52 _ (_ mi - ni -l,-m 2 - n 2 ) 

-L 2 sin (^) e ^2+n 2 )/N sin ^ m Q n _ 52+(-mi - m + 1, -m 2 - n 2 )|, 



where 



(A.13) 



a 1± ^^=— , a 2± ^^=— . (A.14) 



Note that only vectors couple to the ghosts and scalars do not. 
(c) Boson-boson-vector coupling V35 
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V 3b = 2 i ( 2N f /2 V" e -^{mim 2 +nm 2 +(m 1 + ni )(m 2 +n 2 )}/iV- / m ) g M( n ) 
^ Z — ✓ 



<7 



x |Li sin (^-) e-"''"'^ 1 ^ sin (irm on- Trm/iV) 5i + (-mi - n x ,-m 2 - n 2 - 1) 
+Li sin (^-) e" (rai+,ll)/A ' sin (ttto o n + Trni/JV) 5i_(-mi - n x , -m 2 - n 2 + 1) 
-La sin (^) e - 7ri ( ro3+ " 3 ^ JV sin (imon + nn 2 /N) a 2 -(-m x - m - 1, -m 2 - n 2 ) 
_L 2 S i n e ^(™2+« 2 )/JV gin ^ m Q n _ m<2 j N ^ a 2+ (-m x - n x + 1, -m 2 - n 2 )|. 

(A.15) 

(d) four-boson coupling V45 

y « = JT2 E ^(fc)a,(0^(^)«,(«)^ +i+m+ n,oe-" (fclfc2+/li2+mi " l2+ni " 2) 

X ^ (^°^+^ om + mon + no ' e ) _j_ g~ (fcoZ+Zom+mon+nofc) 

(^°^+^ on + nom + mo ^) y (fco/-|-lon+nom+mo/L:) I 



B Cancellation of free energy on torus background due to 
SUSY at 0(JF°) 

First, we consider the cancellation of the free energy on flat fuzzy space backgrounds. On this back- 
ground, the free energy should be zero at least formally. ( The meaning of the word "formally" is 
explained shortly. ) 

For example, the sum of the two-loop diagrams is proportional to 

f 1 2k ■ I 1 

dkdl {k^ + k>.p.(k + iy \- 

Naively this is zero, because 2k ■ I — (k + I) 2 — k 2 — I 2 . However, if the size N of the matrices is 
finite, there should be some momentum cutoff, and (|B.1(I is not necessarily zero. For example, in the 
region where k and I are allowed but k + I is outside the cutoff, the second term of IjB.lfl does not 
make sense. However, such problem occurs only in the large-momentum region. On the other hand 
in the small-momentum region, the free energy vanishes due to the supersymmetry. Therefore, it is 
natural to expect that the free energy is "formally" zero at any order of the loop expansion, if we use 
2k ■ I = (k + I) 2 — k 2 — I 2 and shift the integration momentum not taking into account the cutoff. We 
assume this property. 

Next, we consider the fuzzy torus background. Note that in this case the cutoff procedure is 
completely determined and there is no ambiguity. To be concrete, we consider the T 2 -background 
here. If we use the approximations such as a M (mi, m 2 ± 1) ~ a M (mi, m 2 ), we have 



Via^m) = WNfi(rn)ap(m), (B.2) 
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where 



/i(m) = Li sin 2 (^), 

, . . . /7TTOi\ /7rmi\ 

/ 2 (m) = -I^sin^—jcos^— ], 

/ 3 (m) = -L 2 sin [—) , 

fi(m) = -L 2 sin(^^) cos^^ 2 -) . (B.3) 



Under this approximation, at 0(F°), the only difference from fuzzy plane is that the momentum ^ ± 
is replaced with fj (m) . 

The fj (m) 's satisfy the following set of equations : 

f 2 (m) = £^ sin2 (^)> (B.4) 

i 

2/(m)-/(-n) = / 2 (m) + / 2 (n)-/ 2 (m + n), (B.5) 

which is an analogue of — 2fc • I = k 2 + I 2 — (k + 1) 2 . Using this, we can express the free energy only by 
/ 2 (m)'s. Since the momentum space is completely periodic in the fuzzy torus background, we can shift 
the integration variables without any obstruction, and find that the free energy is zero at this level. 
If we evaluate the free energy without using the approximations, there appear O (-^-corrections. 

C Validity of jF-expansion 

In this section, we discuss the validity of the ^"-expansion. First we consider the case of T 4 . .F's arc 
expressed in the momentum basis such as 

mi —1,7712 — 1)7713 ,771 4 \m\ ,777-2 ,"73 , m 4 

= _ 8i sin eVi/n^Lle-^+m,)/* sin ( ^(m^ m 2 ) 

. (ir{mi-m 2 )\ 

~ Sm { N[ J' (CJ) 

and V 2 is expressed as 

V / 777 1,7772, 7773, 7774; Til, 772, 773 >™4 



i=l,2 V 1 J i=3,4 



; sin 




(C.2) 



Hence, in the generic region, where ^ = 0(1) etc, one ^-insertion amounts to multiplying 
1 < -J— 

N!+N 2 ~ x/iV' 

However, if the arguments in sin are small, the situation changes. For example, if we take Ni <~ 
N 2 <~ >/~N and consider the region where all arguments in sin are of 0(N- a ) (0 < a < ±), then 
4^ - N a -^ 2 . This shows that in the IR region convergence of the ^-expansion becomes worse. 



In the case that Ni ~ N 2 ~ yiV, the free energy does not have IR divergence. Therefore the 
generic region dominates, and the convergence of this expansion is good. Indeed, from the numerical 
calculations of the 2-loop free energy, we can see that one insertion of -^T decreases the power of N by 
|. On the other hand, in the case of T 2 background, the free energy has IR divergence. Therefore IR 
region dominates and the insertion of T does not decrease the power of N. In this case ^"-expansion 
is not reliable. 
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D Tree-level tachyons and extra zero modes 

On fuzzy torus, there exist zero and tachyonic modes, which we list below. 



D.l fermion zero modes 

To be concrete we discuss the T 2 -case here. We use the following 4-dimensional gamma matrices: 

where a l (i = 1, 2, 3) are the Pauli matrices and er 4 = i ■ 1 2 . 

For the left-handed spinor, the Dirac operator is given by 

a 1 + ia 2 , m „ . a 1 — ia 2 , m m , 
<r^ = (Vi-iV 2 ) + (Vi+iPa) 

ct 3 + «cr 4 ,^ ^ , cr 3 - iff 4 _ ^ , 
+ 2 ~ lVi) + 2 ( 3 + 4) 

- 2VN {l 1 [u, •] -l 2 [v-\ •] )■ (D - 2) 

Hence, there are two left-handed zero modes, 

L2V- 1 \ ( L\ 



LiU j and uy y (D.3) 

Similarly for the right-handed spinors we have two zero modes given by 

-L,U ) and (, L 2 V~i 
These zero-modes are not doublers, because their momenta become zero in the large- N limit. 



Ll U J and L 2 V-i ) • 



D.2 Zero and tachyonic modes of the vector 

Although the numbers of zero and tachyonic modes depend on N, they become constant when N is 
sufficiently large. If the ratio r = = 1, there are 4 zero modes and 5 tachyons. On the other hand, 
if r =/= 1, there are no zero modes and 7 tachyons. 

A direct calculation shows that if r = 1 , (Pi, P 2 , — P3, — P4) is an eigenstate with the smallest 
eigenvalue — 16Ai 2 sin 2 (-j^). This mode corresponds to the deformation of the ratio r. 

E 2-loop free energy : 0(JF°) and 0(JF 2 ) contributions 

Here we evaluate the 0{T°) and 0(T 2 ) contributions to the 2-loop free energy on the fuzzy T 4 . For 
brevity of notation, we introduce the following notations 

Pi(m)=-sin (— ] , F 2 (m) = --sin(— Jcos(— J, 

PaM - -1 sin 2 , P 4 (m) = I sin ^ - 



A \ N J ' A \ N J V N 

P 5 (m) = Asm (-jy-J. *e(m) = -Asm (^-^- J cos (^-^-J , 

P 7 (m) = -Asm 2 J , P 8 (m) = Asm J cos J , 

where A = a/ r\fR. We further introduce the abbreviation P(m) • P(n) = ^ • Fj(m)Fj(n), which is 
the analogue of the inner product of momenta. 
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E.l 0(f°) 

We have four kinds of 1PI 2-loop vacuum diagrams ( Fig^l. We list the contributions from each 
diagram : 

(a) Diagram with two fermion-fermion-boson vertices Vzf 

- 2 ^ F(m)-F(-n) 



4N 2 L 4 ^ F(m) 2 ■ F(n) 2 ■ F(m + n) 2 



g 2 i 



E p> (m V.F(r,V (l-cos(27rmon)), (E.l) 



8N 2 L* ^ F(m) 2 ■ F(n) 2 
Here we have used 2F(m) ■ F(—n) = F{m) 2 + F(n) 2 + F(m + n) 2 . 
(b) Diagram with two ghost-ghost-vector vertices V^ g 

The exact result is given by 

* r L 2 sin(^)sin(^) 
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EE 



5127V 2 L 4 ^— ' 1 F(mi + ni ± 1, m 2 + n 2 , + n3, m4 + n4) 2 • F(m) 2 ■ F(n) 2 

: m.fi K 

. . n 2 — m 2 \ ( m 2 + n 2 

x { cos 2nm on± n — cos tt 



Ni J \ Ni 

Lfsin(^) sin(^) 



F(mi + ni,m2 + n 2 ± 1, ^3 + n3, + m) 2 ■ F{m) 2 ■ F(n) 2 

( ( rii-TOi\ / mi+ni 

x < cos 2nm o n =F 7r — — cos n 



Ni ) \ Ni 

+ (L ll N 1 ,mi,m 2 ,n ll n 2 ^ L 2 , N 2 ,m 5 ,m il n d , 1 n i ) . (E.2) 



Shifting m.;'s and introducing approximations, 

1 1 



F(mi + nx ± 1, m 2 + n 2 , m 3 + n 3 , TO4 + 714)2 F(mi + m, m 2 + n 2 , m 3 + n 3 , m 4 + n 4 ) 2 
and so on, we obtain 

g 2 ^—^ 1 — cos(27rm o n) 



(E.3) 



5127V2L 4 ^ F(m) 2 • F(n) 2 
(c) Diagram with two boson-boson-vector vertices Vsb 
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The exact result is given by 



EE 



512iV 2 ^ t—^ F(mi + ?ii =p 1, m 2 + «2, m-3 + "3, m-4 + «-4) 2 ■ F(m) 2 ■ F(n) 2 

( 2 ( T^mi \ f / 27rri2 
x ^ sin si— cos Z7TTO o n T 

. { irm 2 \ . f Tm 2 \ ( f^{m 2 +n 2 )\ ( ir(m 2 - n 2 ) 

sm {-W ) sm llvT J \ cos — ) - cos r m n ± ~~ ^ — 

9g 2 L 2 ^^ 1 



EE 



512iV 2 ^ ^ F{m\ + m, m 2 + n 2 ± 1, m 3 + n 3 , m 4 + ra 4 ) 2 ■ F{m) 2 ■ F(n) 2 

. 9 / 7rm! \ f / 27m! 

x <^ sin <1- cos 2itm o n =F 



AW I V ^1 

. /7rmi\ . /7TOi\ f /V(mi + ni)\ /„ . Tr(mi-ni)' 
— sm sm < cos — cos 2irm o n± 

UJ Ui/l I / v m 

+ (L 1 ,N 1 ,mi,m 2 ,n ll n 2 ^ L 2 , N 2 , m 3 , m 4 , n 3 , n 4 ) . (E.5) 



Again by shifting m and using the approximations (|E.3J) . we can rewrite it as 

27<? 2 1 — cos(27rm o n) 



F(m) 2 -F(n) 2 ' 1 j 

(d) Diagram with two 4-boson vertex I4& 

45g 2 1 - cos(27rm o n) 

256N 2 ^ F(m) 2 -F(n) 2 ' 1 ' ' 

Summing up (a)~(d), we find that the 2-loop vacuum diagrams almost cancel. Actually, if we use 
the approximate values (|E.4|I and (|E .f>|) . we have 



. . . . . . . . / 1 1 27 45 \ a 2 v^l - cos(27rm on) , 

' ' ' r ' ^ ' ^ V 8 512 512 256 / N 2 ^ F(m) 2 ■ F(n) 2 

This cancellation is not exact, but we can see numerically that contribution of the 0(!F Q ) terms is 
smaller than that of the 0{T 2 ) terms. 
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E.2 0{F 2 ) 

The calculation of the 0(J- 2 ) contribution is very tedious. Here we show only the simplest case, that 
is, the planar part of FigQ] (b) : 



g 2 n 2 f 1 - COS { — U ] COS 



E 



ni.i) 



( 2-K{m2+n 2 ) \ 



256L 4 N 3 ^{ F{m) 2 ■ F{n) 2 ■ (F(m + n) 2 ) 3 



1 /. / nmi\ . / nni\ ( n(mi + n 4 )\ . / 7rm 2 \ . f nn 2 \ /n(m 2 +n 2 ) 
sin sin cos + sin sin cos 



1 -cos (^^3)) cos ( 



F(m) 2 • F(n) 2 • (^(TO + n) 2 ) 3 

zr^( ■ { 7rm 3\ ■ ( 7Tn 3\ /7r(m 3 +Ti3)\ . /7rm 4 \ . /7m 4 \ / 7r(m 4 + n 4 ) \ \ 1 

xr'WR sm — — sin — — cos — — + sin — — sin — — cos — — > 

\ \N 2 J \N 2 J V ^2 J \ N 2 J \N 2 ) V ^2 //J 



g 2 ir 2 N f ^ ^ [ 1 - cos (27r(fei + h)) cos (2tt(/c 2 + i 2 )) 



256L 4 7 \ F{k) 2 ■ F(l) 2 ■ (F(k + l) 2 ) 3 

x ^" (sin(7rfci) sin(7rZi) cos (7r(fci + h)) + sm(irk 2 ) sin(7rZ 2 ) cos (ir(k 2 + l 2 ))) 

r 3 yR 

1 - cos (2n(k 3 + l 3 j) cos (2n(k 4 + l 4 j) 
+ F(k) 2 -F(l) 2 ■ (F(k + l) 2 ) 3 

xr 3 VR (sin(7r/c3) sin(7rZ3) cos (7r(fc 3 + l 3 )) + sin(7rfc 4 ) sin(7rZ 4 ) cos ("7r(fc 4 + Z 4 ))) 



(E.9) 



Here, ki — m,-/jVi, kj = rrij/N 2l li — m/Ni, lj = nj/N 2 (i = 1, 2; j — 3, 4), and fc M and vary from 
to 1. 
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